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ABSTRACT 

CTA 102, classified as a flat spectrum radio quasar at z= 1.037, produced 
exceptionally bright optical flare in 2012 September. Following Fermi-L AT de- 
tection of enhanced 7-ray activity, we densely monitored this source in the optical 
and near-infrared bands for the subsequent ten nights using twelve telescopes in 
Japan and South-Africa. On MJD 56197 (2012 September 27, 4-5 days after 
the peak of bright 7-ray flare), polarized flux showed a transient increase, while 
total flux and polarization angle remained almost constant during the "orphan 
polarized-flux flare" . We also detected an intra-night and prominent flare on MJD 
56202. The total and polarized fluxes showed quite similar temporal variations, 
but PA again remained constant during the flare. Interestingly, the polarization 
angles during the two flares were significantly different from the jet direction. 
Emergence of a new emission component with high polarization degree (PD) up 
to 40% would be responsible for the observed two flares, and such a high PD indi- 
cates a presence of highly ordered magnetic field at the emission site. We discuss 
that the well-ordered magnetic field and even the observed directions of polar- 
ization angle which is grossly perpendicular to the jet are reasonably accounted 
for by transverse shock(s) propagating down the jet. 

Subject headings: Galaxies: jets 
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Introduction 



Blazars are highly variable active galactic nuclei (AGN) emitting radiation at all 
wavelengths from radio to 7-rays. They have strong relativistic jets aligned with the 
observer's line of sight and are apparently bright due to relativistic beaming. Their emission 
typically consists of two spectral components. One is attributed to synchrotron radiation at 
lower energies peaking in the radio through optical bands, and the other is inverse Compton 
scattering peaking in the 7-ray bands. Outstanding characteristics of blazars are their 
rapid and high-amplitude intensity variations or flares. These variabilities are observed in 
various wavelengths and time scales. Micro variability (intra-day variability) of flux in the 



optical band has been also detected on timescales as short as minutes to hours (e.g. 



Racine 



1970 



Miller & Noble 



19961 ). It is important to measure the time scales of micro variability 



because it provides limits on the size and location of the emitting regions. 

Polarized radiation is one of the evidences of synchrotron origin in low energies and it 
also varies drastically. Therefore, optical polarimetric observations also pr ovide a strong 



tool to probe jet structures (e.g. 



Marscher et al. 



2008 



Abdo et al. 



2010ah . Nevertheless, 



simultaneous short-term (intra-day) observa tions of flux, color anc 



been performed on 



Sasada et al. 



2008 



y in a few blazars (e.g., 



Andruchow et al. 



Hagen-Thorn et al. 



2008 ; 



Andruchow et al. 



2003 



polarization have 



Cellone et al. 



2007 ; 



20111 ). and hence the origin 



of micro variability is still unclear. 



CTA 102 (al so known as PKS J2232+ 1143, R.A. = 22 /l 32 m 36.4 s , decl. = +11° 43' 50".8 



2000, z=1.037, 



Hewitt &: Burbidge 



( Harris fc Roberts 



19891 ) was first identified as a strong radi o source 



19601 ) and classified as a flat spectrum radio quasar (FSRQ, 



Abdo et al. 



2010bJ) from multi-wavelength o bservations. CTA 102 showed micro variability of optical 



flux and color in the 2004 flare ( jOsterman Meyer et al. 



20091 ). In this flare, the "redder 



when brighter" trend was reported and it was able to be explained by the superposition of 
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variable synchrotron component from radio to optical bands and non-variable "blue bump" 
component which is thought to be thermal disk radiation connecting to the UV band. 



Recently, CTA 



on September 2012 ( iLarionov et al. 



02 showed an extreme activity in the optica 



2012 ; 



Orienti fc D'Ammando 



and GeV 7-ray bands 



2012J). In this paper, we 



present results of high-temporal-density monitoring observations of CTA 102 just after the 
September 2012 7-ray flare. 

2. Observation 

Observations were carried out as a ToO program of Optical and Infrared Synergetic 
Telescopes for Education and Research (OISTER). OISTER is a global observing network 
that links organically many ground-based small telescopes in Japan, South Africa, and Chile 
under a Japanese interuniversity cooperation regime. OISTER is aimed at investigating 
potential transient sources (7-ray bursts, AGNs, supernovae, cataclysmic variables, and 
so on). The biggest advantage of OISTER is its capability of performing continuous and 
high-temporal-density monitoring in many bands extending to as long a wavelength as K s . 
ToO observation of CTA 102 with OISTER was conducted from September 23 to October 
3 in 2012, following the bright GeV 7-ray flare. We obtained the B, V, g\ r\ Rq, Iq, %\ 
z\ J, H, K' and K s band photometric and the Re band polarimetry data with OISTER 
and also with other collaborative telescopes. Telescopes and instruments used for this 
observation are listed on Table [TJ Note that we treated the r' and K' band data as the Rq 
and K s band data, respectively. Reductions of optical and near-infrared (NIR) data were 
performed under the standard procedure of CCD photometry. The position of comparison 
star is R.A. = 22 h 32 m 41.5 s , decl. = +11° 43' 14". 1 (J2000). The magnitudes in the optical 
bands were obtained differentially with nearby compari son stars which have been calibrated 



with the photometric standard stars in Landolt fields (JLandolt 



19921 ) observed on clear 



-7- 



Table 1: Observatories, Telescopes and Instruments 



Observatory/Telescope 



diameter 3- Instrument 



Filters 



Opt/NIR in Japan 



Nayoro Observatory/Pirka 

Akeno Observatory/MITSuME^ 2 ) 

Kyoto University /- 

Koyama Astronomical Observatory/ Araki 

Nishi-Harima Astronomical Observatory/Nayuta 

Bisei Spaceguard Center /- 

Okayama Astrophysical Observatory/- 

Okayama Astrophysical Observatory/MITSuME' 2 ) 

Higashi-Hiroshima Observatory/Kanata 

Iriki Observatory/- 



160 cm 


MSIW 


V,Rc 


50 cm 




g', Rc, Ic 


40 cm 




Re 


130 cm 


ADLER 


B, g>, V, I c , i', z' 


200 cm 


NIC 


K s 


100 cm 


Volante 


r' 


188 cm 


ISLE( 3 > 4 ) 


J, H, K s 


50 cm 




9', Rc, Ic 


150 cm 


HOWPob 5 ) 


V, R c , Bo-Pol 


100 cm 


Infrared Camera 


J, H, K' 



Opt/NIR in South- Africa 


South African Astronomical Observatory/IRSF 140 cm 


SIRIUS 


J, H, K s 


Radio Observatory 


Mizusawa VLBI Observatory/Hitachi 32-m Telescope 32 m 




8.4 GHz 



a Size of primary mirror 



References. 



;D 



Watanabe et al 



Yanagisawa et al. 



20081 : (5) 



Kawabata et al. 



2012; (2) 



Kotani et al 



2005 



2008. 



(3) 



Yanagisawa et al. 



2006; (4) 



and stable nights. The same comparison st ars are used also in N IR bands, where the 

i). We corrected the data 



magnitude is given in the 2MASS catalogue 
for the Galactic extinction (e.g., ^=0.233, 



fSkruts kie et al. 



2006 



Schlafly fa Finkbeiner 



2011 



NED database^!). 



There were small systematic difference in photometric system among observatories and 
instruments; the standard deviations of the magnitudes of the comparison star during the 
observation period were ARq ~ 0.02 mag (~ 2% of flux) and AK S ~ 0.06 mag (~ 5% of 
flux). These values were added to the photometric errors of CTA 102 in each band. 

The polarimetric observations were performed wi th HO WPol installed to the Kanata 



2008|). A unit of the 



Telescope located on Higashi-Hiroshima Observatory (IKawabata et al. 
observing sequence consisted of successive exposures at 4 position angles of a half-wave 
plate; 0°, 45°, 22.5° and 67.5°. Polarimetry with HOWPol suffers from large instrumental 
polarization (Ap ~ 4%) produced by the reflection of the incident light on the tertiary 
mirror of the telescope. The instrumental polarization was modeled as a function of the 
declination of the object and the hour angle at the observation, and we subtracted it from 
the observation. We estimated that the error in this instrumental polarization correction is 
smaller than 0.5% from many observations for unpolarized standard stars. The polarization 
angle (PA) is defined in the standard manner as measured from n orth to east. The PA 



19831 ). 



was calibrated with two polarized stars, HD183143 and HD204827 (ISchulz fa Lenzen 
Because the PA has an ambiguity of ±180° x n (where n is an integer), we selected n which 
gives the least angle difference from the previous data, assuming that the PA would change 
smoothly. The error of PA was estimated to be smaller than 2° from observations of the 
polarized stars. 

The radio data were obtained from September 26 to 28 using Hitachi 32 m telescope 
of Mizusawa VLBI Observatory, NAOJ, which is operated by Ibaraki University. The front 



1 http://ned. ipac.caltech.edu/ 
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end was a cooled HEMT receiver, and a typical system temperature was 25 K including 
the atmosphere toward the zenith. Since beam switching system is not equipped, we 
rapidly scanned the antenna around the target CTA 102 in azimuth and elevation direction, 
while recording the total power. As a result, fluctuation of the observed power due to the 
atmosphere can be minimized and the pointing error can disappear. The accuracy of the 
calibration was estimated to be 10%. 

3. Results 

Figure [1] shows temporal variations of the optical i?c _ band and NIR i^ s -band total 
fluxes, as well as those of the i?c-band PD and PA during our dense monitoring of CTA 102, 
following intense 7-ray flare detected by Fermi-L AT. In the third panel of Figure 1, we also 
show a light curve of i?c _ band polarized flux (PF), which was calculated by 

100 K J 

where PD and Fr c are the measured polarization degree in unit of % and total flux in the 
Re band. Except for some nights with bad weather, we were able to monitor CTA 102 
continuously. We note that, although not shown here, light curves in other bands exhibited 
quite similar temporal variations to those in the Rq and K s bands. 

First, the light curves in the Rq and K s bands showed a clear decay on MJD 56194. 
Given that Fermi-LAT detected a strong 7-ray flare on MJD 56189, the observed declining 
profile likely corresponds to the decay phase of the bright 7-ray flare. After that, both of 
the light curves showed a mild and symmetrical enhancement peaking around MJD 56198 
until the end of this follow-up observation, except an intra-night strong flare on MJD 
56202. In addition, we can clearly see some remarkable features in the daily polarimetric 
data. Namely, possible PA swing from ~ 0° to ~ 100° was observed on MJD 56195 and 
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56196. To check the temporal evolution of the PA, we plotted the Stokes parameters Q and 
U observed during these nights in Figure |21 together with those of other nights. Again, 
gradual rotation of PA was obviously seen, confirming the PA swing by ~ 100° during these 
two nights. 

The continuous PA rotation was terminated by a sudden jump to ~ 50° on MJD 56197. 
More strikingly, the PF rapidly increased by three times from ~ 0.1 x 10~ n erg cm' 2 s _1 to 
~ 0.3 x 10 -11 erg cm -2 s _1 without an apparent increase of the total flux during the night. 
This is the first clear detection of a short-term (hours-scale) PF flare without corresponding 
total flux enhancement. Another interesting feature is a strong intra- night flare detected on 
MJD 56202. The isolated flare showed smooth and symmetrical profile peaking around 13 
UT, and the total and polarized fluxes increased by factors of ~1.5 and ~2.0, respectively. 
Notably, the PD evolved in a quite similar way as the total flux, while the PA remained 



almost steady during the flare (see also Figure [Hand Figure [2]). The PA m easured on M 



56202 is orthogonal to the jet direction observed from VLBA observation (JFromm et al. 



D 



20121 ). This is in contrast with the case of AO 0235+164, in which the 



optical PA aligned 



near 



2008|). 



y the same direction as the jet during the flare in 2006 December (IHagen-Thorn et al. 



To see spectral evolution in the optical and NIR b ands, we plotted Re, — K K colo r 



Qsterman Meyer et al. 



( 120091 ) . a 



against i?c-band magnitude in Figure [3j As reported in 
"redder when brighter" trend is excepted in micro variability. As a result, however, there 
is no clear correlation between Re — K a color and i?c-band magnitude for a whole data. 
Looking at the colors measured during our observation other than MJD 56197 and 56202, 
the source showed a hint of "redder when brighter" trend, as is reported previously by 



Qsterman Meyer et al 



(120091) 



On the other hand, radio fluxes were almost constant within 10% during our monitoring 
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and were comparable to that measured in the quiescent state (iStanghellini et al.l 119981 ). 
This is because the emission region would be optically thick at 8.4 GHz. Actually, the 
const ant radio flux durin g optical and 7-ray flare was previously observed for CTA 102 in 



2006 (IFromm et al. 



201 lh . 



4. Discussion 



On MJD 56197, a transient PF flare without significant increase of the total flux was 
observed. This "orphan PF flare" is interesting and explained in such way that a new 
emission component which is less luminous but possesses extremely high PD with respect 
to the long-term baseline component might have emerged suddenly. Large discrepancy of 
PA between MJD 56196 and 56197 (~ 100° to ~ 45°) supports the emergence of the new 
component. A careful inspection of the light curve showed that the source brightened only 
by ~10% in the total flux (which corresponds to ~ 0.2 x 10~ n erg cm -2 s _1 ) during the 
intra-night increase of PF from 0.1 x 10~ n erg cm -2 s _1 to 0.3 x 10~ n erg cm -2 s _1 , implying 
that the new component has the PD of ~100%. To overcome this unrealistic situation, 
we have to consider that decrease of total flux of the long-term component occurred 
simultaneously with the emergence of the highly-polarized new emission component and 
weakened increase of the total flux. In this case, the PD of the new component depends 
on how much total flux of the long-term component decreases. For example, we assume 
that gradual decrease of total flux of the long-term component by 0.2 x 10~ n erg cm -2 
s _1 occurred simultaneously with gradual increase of the new component by 0.4 x 10~ n 
erg cm -2 s _1 Then, the PD of the new component can be estimated as 50%, which seems 
reasonable compared to the above ~100% PD scenario, though still very high. We note 
that some blazars such as O J 287 and BL La c exhibited similar orphan PF flares in 



longer day-month timescale (jlkejiri et al 



20111 ). It is currently unclear whether the same 
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Fig. 1.— Multi- wavelength light curves of CTA 102 from September 24 to Octobe r 3 in 
2012. Top panel: NIR (A^-band) flux. Second panel: optical (.Rc-band) flux. Third panel: 
polarized flux (PF) in the Rq band. Fourth panel: polarization degree (PD) in the Rq band. 
Bottom panel: polarization angle in the Re band. 
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mechanism is responsible for such a similar phenomenon with longer timescale. Anyway, 
our observation presented here is the first measurement of the orphan PF flare in much 
shorter (intra-night) timescale, and further theoretical study is needed. 

Smooth and continuous change of PA from ~ 0° to ~ 100° observed o n MJD 56195 



and 5 6196 suggests either the pr esence of he 



20081 ) or global bending of a jet (jAbdo et al 



ical ma gnetic field inside a jet ( iMarscher et al. 



2010al ). Note that the PD remained constant 



at ~ 3% during the PA swing, which was also observed in both of the two scenarios from 
BL Lac (helical magnetic field) and 3C 279 (bent jet). The observed large PA swing rate of 
~ 50° per day is similar to the case of BL Lac and four times larger than that of 3C 279, 
implying that the helical magnetic field scenario might be preferred. On the other hand, if 
PA swing in the opposite direction is observed, the bent jet scenario would be validated. 

Current Very Long Baseline Interferometry (VLBI) polarization observation is now 
feasible to examine magnetic field structure of the parsec-scale jet. Especially, measurement 
of Faraday rotation is recognized as a powerful tool to manifest the presence of helical 
magnetic field. When polarized radiation propagates within a magnetized plasma, the plane 
of linear polarization rotates and the rotation measure (RM, rotation of the polarization 



angle) depends on the line-of-sight component of the magnetic field (see e.g. 



Chen 



19741 ). 



Hence, if helic al magnetic field is surrounding the jet, RM gr adient across the jet s hould be 



observed (e.g. 



Blandford 



19931 ). Indeed, several authors (e.g. 



Gabuzda et al 



20041 ) claimed 



detection of RM gradient transverse t o the VLBI jet in various blazars since the first 
discovery from 3C 2 73 by 



is still controversial ( Taylor fc Zavala 



Asada et a 



(2002), although the significance of these detections 
2010)). Importantly, recent Very Long Baseline Array 



(VLB A) observation detected significant RM gradient from CTA 102 (IHovatta et al.ll2012l ). 
The presence of helical magnetic field is consistent with our observation of relatively fast 
PA swing over the two nights on MJD 56195 and 56196. Thus, optical polarimetry provides 
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another powerful tool to claim the presence of helical magnetic field. 

A prominent brightening was observed on MJD 56202, as is evident from the light 
curves of both the total flux and PF, with almost constant PA of ~ 70° (see Figure 1). 
This could be again explained by emergence of a new bright component with high PD 
with respect to the long-term baseline emission component. Comparison of the total and 
polarized fluxes between the beginning and peak of the flare allows us to estimate the PD 
of the new component. Namely, ratio of the PF increase by ~ 0.35 x 10 -11 erg cm -2 s _1 
with respect to the total-flux increase by ~ 1.0 x 10 -11 erg cm -2 s _1 leads to the PD of 
~30-40% for the new emission component. Such a high PD indicates presence of highly 
ordered magnetic field at the emission site, which might be generated through compression 
of turbulent magnetic fie ld by shocks inside a jet as advocated by previous papers (e.g., 



Hagen-Thorn et al. 



2008|). 



The optical PAs measured during the two flares on MJD 56197 and 56202 were 
45° and ~ 70°, respectively. Compared to the VLBA jet image (see e.g. Fig. A. 15 of 



Fromm et al 



2013 



for the jet image), we found that they are nearly orthogonal to the 
jet direction. Since magnetic field direction is in principle assumed to be perpendicular to 
PA, the measured PAs result in a claim of magnetic field orientation aligned with the jet. 
How is the highly-ordered magnetic field aligned with the jet generated? At first glance, 
the shock-in-jet scenario considered above seems unreasonable, because the magnetic field 
compressed by shocks propagating down the jet is aligned with a direction transverse to the 
jet. It should be noted that here magnetic field direction is inferred based on the assumption 
that it is orthogonal to the observed PA but recent theoretical work sugges ts that such an 



assum ption is not co rrect for relatiy i stical 



y moving sources like AGN jets (jLyutikov et al. 



20051 ). In particular, ILyutikov et al.l ( 120051 ) pointed out that toroidally dominated magnetic 



field is observed as poloidal in the observer's (rest) frame when ultra-relativistic jet with 
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T >> 1 is viewed at small angle to the line of sight (T is bulk Lorentz factor of the jet). 
Given this complicated situation, the measured EVPAs significantly different from the jet 
direction can still be accounted for by the "shock-in-jet" scenario. 

It is interesting to compare the PAs observed during the t wo flares on MJD 56 1 97 an d 



Hagen-Thorn et al. 



(120081 ) 



56202 with past observation, although only one observation by 
is found in the literature. The authors reported hours-scale short timescale polarimetric 
variability from BL Lac object AO 0235+164 during an outburst in 2006 December and 
fou nd that the PAs tend to a lign with the jet direction around the maximum PD (see Fig. 



4 of 



Hagen-Thorn et al 



120081 ). This is quite different from our observation of perpendicular 
EVPAs to the jet direction. If relativistic effect hypothesis is true, uniform distribution of 
EVPA with respect to the jet direction would be found for hour-scale flares. In any case, 
the current samples are only two and further optical polarimetric observation is definitely 
needed. 

To summarize, we performed dense optical/IR photometric and polarimetric monitoring 
of CTA 102 following strong 7-ray flare in 2012 September with OISTER program. We 
found (i) smooth and gradual PA swing by ~ 100 deg over the two nights on MJD 56195 
and 56196, (ii) orphan polarized flux flare on MJD 56197, (iii) significant brightening 
on MJD 56202 and (iv) grossly perpendicular PAs to the jet direction during the two 
flares. Combined with recent VLBA detection of RM gradient across the jet, we infer 
that helical magnetic field would be present at the emission region responsible for the 
long-term baseline component. The observed two flares can be explained by emergence of 
new emission component which possess highly-ordered magnetic field. Such a magnetic field 
configuration would be generated through compression by shocks propagating down the jet. 
The observed EVPAs perpendicular to the jet direction is not unreasonable given the effect 
of relativistically-moving radiation source. 
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